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ABSTRACT 

During atmospheric entry of planetary probes, the thermal protection system (TPS) of the probe 
is exposed to high temperatures under low pressures. In these conditions, carbonous fibrous 
TPS materials may undergo oxidation leading to mass loss and wall recession called ablation. 
This work aims to improve the understanding of material/environment interactions through a 
study of the coupling between oxygen transport in the Knudsen regime, heterogeneous 
oxidation of carbon, and surface recession. A 3D Random Walk Monte Carlo simulation tool is 
used for this study. The fibrous architecture of a model material, consisting of high porosity 
random array of carbon fibers, is numerically represented on a 3D Cartesian grid. Mass 
transport in the Knudsen regime from the boundary layer to the surface, and inside this porous 
material is simulated by random walk. A reaction probability is used to simulate the 
heterogeneous oxidation reaction. The surface recession of the fibers is followed by front 
tracking using a simplified marching cube approach. The output data of the simulations are 
ablation velocity and dynamic evolution of the material porosity. A parametric study is carried 
out to analyze the material behavior as a function of Knudsen number for the porous media 
(length of the mean free path compared to the mean pore diameter) and the intrinsic reactivity 
of the carbon fibers. The model is applied to Stardust mission reentry conditions and explains 
the unexpected behavior of the TPS material that underwent mass loss in volume. 
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. Scientific context 


Pyrolysis and ablation of PIC 


aB 


(Phenolic Impregnated Carbon Ablator) 



Current models 


I New approach 
I for porous media 





SEM micrographs W 
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. Literature data 

Stardust Post-flight Analyses W 
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-> Occurrence of a volume phenomenon in the char layer? 

- Oxidation (oxygen from the atmosphere) 

- Sublimation 

- Reactive pyrolysis gases (H 2 0, 0 2 , C0 2 , ...) 

- Mechanical erosion of the matrix 
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. Ablation model : multiscale approach 

Microscopic simulation of the oxidation of a 


carbonized PICA 


• Hypotheses : 

- T = 3360 K (isothermal) 

- P = 0.26 Atm 

- Air 

- No pyrolysis gas blowing 

- p f = 32 p m = 1 .8 g/cm 3 
[Stackpoole et ai, 2008] 

- k m = 10 k f = 13.7 m/s 
[Drawin 1992, Lachaud 
2007] 

- Simulation : 1.2 s 
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Carbon preform density (kg/m 3 ) 


. Application : Stardust Reentry Conditions 

End of the ablative part of Stardust trajectory : 



90s CW 1350 K < T sublim ) to 130s (T M =850 K <T oxi ) 


Pure oxidation of the char layer of a PICA 
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Loosely coupled with 
traditional tools: 

• data from FIAT inside 
the material 

- Temperature gradient 

- Pyrolysis gas velocity 
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. Application : Stardust Reentry Conditions 

Fit of experimental data using the 



volume ablation model 
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<J) Stardust Post-Flight Analyses 
[Stackpoole, 2008 ] 
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(this study) 
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. Conclusion 



• Multi-scale model for the Ablation layer 

- Surface / Volume ? -> depends on conditions 

- Transition criterion : Thiele Number » 1 -> Surface 

- Stardust reentry conditions: Volume oxidation under the hypotheses of the model 

- First material response model accounting for volume ablation 


-> capability to reproduce experimental data 



“Emerging part of the Iceberg” 


October 1 5, 2008 


Hypersonics Materials and Structures - Ablators Discussion 


7 



. Perspectives 

Integration of volume ablation in the pyrolysis model 
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. Modeling conduction in porous media 

Illustration 
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. Determination of the effective diffusion coefficient 

Direct Simulation Monte Carlo (DSMC) 



Monte Carlo Simulation : 

- Random Walk rules : 

• (T,P) fixed = (A, D) fixed 

• A : Maxwell-Boltzmann distribution 

• constant velocity norm (D = 1/3 v A) 

• 3D random direction drawing 


1) Displacement ^of 10000 walkers 
followed during x (chosen for convergence) 

2) Einstein relation on diffusion process : 



• Tortuosity as a function of Kn for the fibrous preform of this study 



Illustration : path of a walker 
in a periodic cell 
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. Modeling advection in porous media 

Traditional CFD method - 



DSMC to be used for Knudsen regime 



SEM observation Fibrous media geometry model 



€ g v g = 


K, 




Vp, 


K g (virgin)= 1.6 10 ''m 2 


K g (pyrolysed)= 2.0 10' 11 m 2 


Same order of magnitude as the 
Kozeny-Carman law 
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. Perspectives : Summary 
Multiscale modeling strategy 
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Macroscopic model : Ablation & Pyrolysis 

• 3D Cartesian code (Finite Volume / 3D Level set) 

• Experiments 


Microscopic approach : multi-physics 
• Analytic, CFD & DSMC 
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. 3D simulation tool : AMAt 1 ) 
Random walk algorithm / 


D Marching Cube front tracking 


M 



Boundary 
condition 
(Buffer zone} 

Mass transfer 
(Brownian motion) 


Ablation rate 
(Sticking event) 


Moving fluid/solid 
interface 
(simplified 
Marching cube) 


(1) J. Lachaud, G.L. Vignoles. A Brownian 
motion technique to simulate gasification 
and its application to C/C composite 
ablation. Computational Materials 
Science, 2008. In press. 
doi:10.1016/j.commatsci.2008.07.015 
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. Ablation model : multiscale approach 

Effective diffusion coefficient 


• Knudsen effects 

Bosanquet model 

_L._L + JL__L + _i_ (harmonic average ) 

D ref D B D K A V ^ 'A vd p 

E.g. Stardust peak heating Kn ~ 0.1, i.e. A » 10 d 



• Tortuosity effects 

Averaging method 



Tortuosity has to be assessed 
by numerical simulation 
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. Ablation model : multiscale approach 

Determination of the effective diffusion coefficient D 


eff 



Monte Carlo Simulation : 

- Random Walk rules : 

• (T,P) fixed = (A, D) fixed 

• A : Maxwell-Boltzmann distribution 

• constant velocity norm 

(D = 1/3 v A) with 3D random 
direction drawing 




1) Displacement ^of 10000 walkers 
followed during x (chosen for convergence) 

2) Einstein relation on diffusion process : 
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Tortuosity as a function of Kn for the fibrous preform of this study 



Illustration : path of a walker 
in a periodic cell 
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3 kinds of Experiments 

Improvement & Validation of the models 


Device 


Measured data 

Dimension 
& scale 

Time scale 

Difficulty 

Interest 

Other 

comments 

<£EM* 

^Scanning 

Electron 

Microscopy) 

4-« 


Architecture 

(intuitive) 

resolution: 100 nm 

2D surface 
micro to 
macro 

1 month 

Image 

analysis 

Help into 
architecture 
modeling at 
all scales 

Relatively low 
cost for fast 
preliminary 
results 
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Chemistry 
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these gas 

ID 
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- No rush 
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1970 

- Useful 
again for 
porous TPS ; 

jSihort ramp 
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0-2500 K 


Effective conductivity 
k=f(T) 

3D 

macro 

6 month 
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e position 
-data analysis 

Useful for 

radiation 
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Needed 

ASAP 

Long steady 
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2+ tangential 
(flowing) 


Thermal gradient 
Density profile 
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3D micro 
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ple position 
-quantify 
spallation= f( 
shear stress) 

Ablation/pyro 
lysis coupling 
Spallation 
Model 
validation 

Parametrical 
study on 
blowing and 
grad(T) 

Plasma 
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idem 

3D micro 
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macro 
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-test 

conditions 
-data analysis 

Global vali- 
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fluid + mater 

Parametrical 
study on P, 
T, v 


* : to be done in priority with available funding / ** : to plan now and begin in 1 year 
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